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Abstract 

We study viscous hydrodynamics of hot conformal field theory plasma with 
multiple/non-Abelian symmetries in the framework of AdS/CFT correspondence, 
using a recently proposed method of directly solving bulk gravity in derivative 
expansion of local plasma parameters. Our motivation is to better describe the 
real QCD plasma produced at RHIC, incorporating its U(l) N f flavor symmetry as 
well as SU(2)j non-Abelian iso-spin symmetry. As concrete examples, we choose 
to study the STU model for multiple U(l) 5 symmetries, which is a sub-sector of 
5D N=4 gauged SUGRA dual to N=4 Super Yang-Mills theory, capturing Cartan 
C/(l) 3 dynamics inside the full R-symmetry. For SU (2), we analyze the minimal 4D 
N=3 gauged SUGRA whose bosonic action is simply an Einstein- Yang-Mills system, 
which corresponds to SU(2) R-symmetry dynamics on M2-branes at a Hyper-Kahler 
cone. By generalizing the bosonic action to arbitrary dimensions and Lie groups, we 
present our analysis and results for any non-Abelian plasma in arbitrary dimensions. 
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1 Introduction and summary 



Strongly coupled plasma of finite temperature gauge theories has recently become a fasci- 
nating subject of research, largely motivated by the RHIC experiment of relativistic heavy 
ion collisions. Naive QCD expectations based on perturbative QCD have failed to explain 
certain important aspects of the created QCD plasma, and there are several indications 
that the RHIC plasma is in fact a strongly coupled liquid. Given the situation, one may 
hope that the problem can be attacked by AdS/CFT correspondence or gauge/gravity 
correspondence because the correspondence is useful precisely when the gauge theory side 
is strongly coupled pQ. In the gravity side, a finite temperature plasma corresponds to a 
black-hole, or more precisely black-brane, spacetime with Hawking temperature identified 
with the temperature of the gauge theory. The black-hole horizon is located at certain 
point in the holographic additional dimension, and presumably physics outside the hori- 
zon with suitable boundary conditions on the horizon describes the finite temperature 
plasma of gauge theories. There have been a lot of useful and often surprising results 
obtained from this gravity picture, which would be hard to be found in the pure gauge 
theory analysis due to strong coupling [2j [3j [U [H [91 [10] . 

Of course, the main huddle that hinders further progress in this direction is the absence 
of precise dual gravity theory of the real QCD at present, although it is important to 
improve the current models to better mimic realistic QCD (TTJ [12J, [131 E]- However, 
certain properties of strongly coupled finite temperature plasma may be universal at least 
qualitatively [I]; a well-known example is the viscosity-entropy ratio [21 HUH], r\js ~ 4-, 
which in fact is close to the RHIC experiment data, and it is not a vividly wrong idea to 
try to learn something about realistic RHIC plasma by studying certain specific AdS/CFT 



models at finite temperature 3. The question whether the results obtained in the specific 
model are meaningful in realistic QCD should be asked carefully though. 

Any finite temperature plasma is described by hydrodynamics in sufficiently slowly- 
varying and long-ranged regime. It is more a framework rather than a result; it is based 
on local thermal equilibrium and conservation of symmetry, such as energy-momentum 
or global symmetry currents. It is then natural to expect that hydrodynamics should 
be emerging in the gravity side of finite temperature plasma described by black-branes. 
Indeed, it is based on this idea that various hydrodynamic coefficients including the 
viscosity-entropy ratio were calculated via linear response theory in the gravity back- 



1 See [5] for possible violations of viscosity-entropy bound. 
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ground [21 [3j HI [6j [7J [TO], [15]. ^46 m#io way of deriving the hydrodynamics from the 
gravity side when black-brane parameters like horizon and charges are slowly varying was 
recently developed (TBI [T71 [HI [19], and extended to the single U{1) R-charged system 
in ref.[20l |2T| [22]. This progress is important because one can in principle go beyond 
the linearized approximation to arbitrary non-linear orders one desires, and some of high 
order transport coefficients and non-linearity seem interesting |25| . For applications to 
non-relativistic AdS/CFT, see ref . [26j [27J [28] , and for dyonic system, see ref . J29J, [30] . See 
also ref. [3TI [32] for similar developments. 

In this work, we generalize this line of development in two different ways, motivated 
by real QCD plasma. In QCD with several quark flavors, there is an enlarged global 
symmetry SU(Nf) L x SU(Nf) R x U{1)b- If one first focuses on the quark species, each 
quark flavor has its own conservation; in the case of three quarks u, d, s that seem relevant 
in the RHIC experiment, one should deal with finite temperature plasma with U(l) 3 global 
symmetry. Note that these U(l) 3 components are highly interacting with each other 
by strong interactions, and except their conservation laws one can not predict a priori 
anything about their dynamics such as diffusion coefficient etc. Because these interactions 
are crucially affecting the results of transport coefficients, we better work in a well-defined 
AdS / CFT set-up rather than working in an arbitrary unguided gravity theory. We choose 
to study the STU model, which is a consistent truncation of AdS$ x S 5 with U(l) 3 (or any 
Toric Sasaki-Einstein compactification) dual to N = 4 SYM plasma with three Cartan 
C/(l)'s inside SO(6)r, as a model example of multi-charged finite temperature plasma. 
This model has been previously studied in linear response approach in ref. [23, 24J, but our 
framework enables one to go beyond linearized approximation. We also compute charge 
diffusion coefficients in the model for the first time in the literature. Although our set-up 
is not precisely QCD, we hope that it captures some aspect of real U(l) 3 dynamics in 
RHIC plasma. 

Our second generalization is for non-Abelian SU(2), although our result is valid for 
an arbitrary Lie groupB This has a clear motivation from QCD again; it mimics iso-spin 
SU(2)v symmetry of QCD in two-flavor approximation in a late stage of RHIC plasma. As 
fluctuations of charged pion density correspond to SU(2) V fluctuations, our study should 
be interesting in describing pion fluid at finite temperature too. For a specific model, we 
study a realization of SU(2) gauged supergravity in Tri-Sasakian compactification of M- 
theory to AdS^, corresponding to the S , L r (2)^-symmetry sector of 3-dimensional M2-brane 

2 For study of SU{2) in a different context of condensed matter system, See ref. [33l l34l [35] . 
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plasma. However as we are more interested in higher dimensions such as 4-dimensions, 
we simply generalize the bosonic action to arbitrary dimensions AdS n+ \ with n > 3 
and perform analysis in complete generality of dimensions. Our analysis automatically 
includes the generalization of the previous single U(l) case to arbitrary dimensions as well. 
One should however note that our U(l) 3 and SU(2) R-symmetry in the field theory side 
do not come from fundamental flavors, but rather from flavors with adjoint representation. 
To introduce fundamental flavors, one normally needs to introduce extra branes in the 
set-up, whose detailed study in hydrodynamics is remained for future work. 

We obtain the results at first order in derivatives, while the second order calculation 
is straightforward, as we will mention in the text, but extremely complicated to present. 
We leave its more controlled analysis to the future. We however mention that even the 
first order transport coefficients we obtain are quite non-trivial. In the STU model, we 
find three conserved currents of each U(l) 3 at first order in derivatives to be 

J? = pm* - Vj (rT + D uPl + de v ^u v d p u a + ■■■ , (1.1) 

where I = 1,2,3 runs for U(l) 3 symmetries, p 1 is the charge density, and the diffusion 
coefficients Vj are given by 

V,= (ri ~ gj) , (1.2) 

with r H being the horizon radius, and the relation of qj with the energy/charge densities 
can be easily found in the text. The parity- violating coefficients Q originated from the 
5D Chern-Simons term are 

j- 1 ( n y/rnqjy/mq K 



32nG 5 \ ^K + qjXrjj + qK) 



VrnqJ r jmq]^mq K jmqi \ 

3m ^^(rl + qj^rl + qK^rl + q,)) ' 1 ^ 

where C IJK is the Chern-Simons coefficient. As far as we know, this is the first time in 
literature to have these results. For the SU(2) case, our result for the first order correction 
to the SU (2) currents in n-dimensional CFT is 

Jf ] = -v ( p ' P ^ dvp) - u»d v u\ P a + v x ^p h p;{d v p c ) 

\ P P / 

+ v 2 p;{p a {p.d u p)-{p-p){d uP a )) , (i.4) 

with three diffusion coefficients 

v = 1 L _ 2q- q \ = (n - 2)m + 2r\ 

(n — 2)r H \ nmr 1 ^ 2 J n(n — 2)mr H 
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87rGWi/ (n ~ 2) 
(n-1) 



Z>2 



3 1 

(n — 1) 2 (n — 2) 2 



(1.5) 



The £> is essentially the usual diffusion coefficient of Abelian nature, while the other two 
diffusion coefficients are due to the non-Abelian properties. The constants /( n ~ 2 ) and 
g( n ~ 2 ) are defined in the text. We hope that this structure is a useful starting point to 
study non-Abelian iso-spin plasma of RHIC. 

2 Crash review of the method 

The basic idea in ref. [HJJ [TTj for deriving hydrodynamics from gravity is conceptually 
quite neat; given a black-brane solution with certain parameters such as temperature, 
charges etc, one simply considers slowly-varying those parameters in the solution. Note 
that keeping the form of the original solution while only varying the parameters inside 
has a physical meaning of local thermal equilibrium with given parameters at that point. 
However, the resulting configuration with varying parameters will no longer solve the 
equations of motion by obvious reason, and one should add corrections to the original 
form of the solution to satisfy the equations of motion. These corrections will clearly 
be sourced by the derivatives of the black-brane parameters one is slowly-varying, and 
one can systematically invoke derivative expansion for these corrections. After obtaining 
the full solution at fc-'th order in derivatives, one can read off physical quantities at that 
order via AdS/CFT dictionary, such as energy-momentum tensor and charge currents. In 
principle one can go to an arbitrary order in derivative expansion systematically. 

We simply illustrate the procedure for the metric and we refer to the original work of 
ref. [T6l [T7] for more complete discussion. Suppose we have a homogeneous black-brane 
solution with the metric gf^v(u M , m, Qi) where w M , m and Q^s are parameters of the 
solution such as 4-velocity, energy density and charge density of the plasma. Once we 
allow for these parameters to vary in CFT coordinate directions up to fc'th derivative, one 
should add correction terms 



to solve the equations of motion, where . . . are varying in CFT coordinate x M up to 

fc'th derivatives only, and g^ N is a local function of derivatives at i'th order. Suppose one 
solved this problem up to fc'th order. Then to go to the next order, one simply considers 
the varying parameters in the above k 'th solution up to (k + l)'th order, which would not 



k 




(2.6) 



i=l 
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solve the equation of motion any more due to (k + l)'th order in derivatives we are now 
considering. Because the equations of motion are solved up to fc'th order already by the 
above, one simply needs to add g^^ to 9mn that is a local function of (k + l)'th order 
in derivatives. Typically the equations of motion are second order partial differential 
equations, and since the variation of itself along would be the next order to be 

neglected, g^M^^ only depends on the holographic direction r and the local derivatives of 
parameters, without any x M dependence at this order. Therefore one would get a simple 
second order ordinary differential equation along r direction from the equations of motion, 

Lr {g^MN ^ = &MN ^ 1 (2-7) 

with a source S^rP being some function of local derivatives at (k + l)'th order. Note 
that L r would be universal, without being dependent on k, completely determined by 
the zero'th order solution, and one can obtain the source S^iP quite straightforwardly 
by plugging the above fc'th order solution with (k + 1) 'th order derivatives of parameters 
into the equation of motion, and simply gathering uncanceled (k + l)'th order term left. 
Therefore even after performing the analysis at the first order, one can find the ordinary 
differential operator L r , and the subsequent higher order analysis will then become con- 
ceptually simpler. One can go on these steps inductively to arbitrary order in derivatives. 



3 Hydrodynamics with U(l) 3 : the STU model 

The action of the STU model which is a sub-sector of AdS§ gauged N = 4 supergravity 
holographically dual to iV = 4 SYM theory is 

(16ttG 5 )£ = R + 2V(X) - l -G IJ {X){F I ) MN {F J ) MN -G IJ {X)d M X I d M X J 

+ ^T^e MN ^ R C I jK{F I ) M N{F J ) PQ {A K ) R , (3.8) 

where 

and Cuk is totally symmetric with C123 = 1. Also, X 1 are not independent but con- 
strained by 

hj IJK x I x J x K = X 1 X 2 X 3 = 1 . (3.10) 
We have put L 4 = 4ng s Nl 4 = 1 for simplicity and in this convention, we have 

G 5 = ^ , (3.11) 
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where N is the rank of the gauge group. Capital letters M, N, ■ ■ • represent 5-dimensional 
indices, while Greek letters fi,v,- ■ ■ would mean 4-dimensional indices. The equations of 
motion one obtains consist of the Einstein equation 

Rmn + (~V(X) + \g ij (X)(F i ) mn (F j ) MN )9mn 

- G IJ {X){F I ) PM {F J ) P N - G /J (X)9 M X / ^X J = , (3.12) 

the three Maxwell equations for each I = 1, 2, 3, 

V N (G LJ (X)(F J ) MN ) - — ^e MNp Q R C IJK (F J ) NP (F K ) QR = (3.13) 
' 16V— (7^ 



and the scalar field equations 



~\ (^W 1 ) + l(F'MF J ) MN ) % = , (3.14) 

where 0* (i = 1, 2) are any independent parametrization of X s. 

The black brane solution with arbitrary three charges has been known [36] and is given 

by 



ds = -H s {r)f{r)u lM u v dx^dx v - 2H e (^u^dx^dr + r H 3 (r) (77^ + w^) dx^dx v 

A- = , X'=?*> , (3.15) 

r z + gj ni[r) 

where 

■? 

/(r) = --+r 2 #(r) , H(r) = l[H I (r) , F 7 (r) = 1 + ^ , (3.16) 
r I=1 r 

and m m is the 4-velocity of the fluid with u^u^ = —1. Our convention is r]^ v = diag(— , +, +, +) 
and = (—1, 0, 0, 0) in the rest frame. 

As we discussed in the previous section, we consider slowly varying parameters u^, 
m, and qi up to first order, and we work in the frame where = (— 1, 0, ■ ■ • , 0) at the 



3 The gauge fields in the solution have an infinite norm at the horizon due to the diverging g 00 , which 
can be remedied by going to the grand-canonical ensemble with chemical potential that can be added to 
the solution as a constant mode. We thank Mukund Rangamani for pointing this issue to us. 
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position x M = for simplicity. Once we find the solution, we can easily make the result 
relativistically covariant. Then at first order in derivatives, we have 



m 
Qi 



+ x^d^m 



(3.17) 



and the above black-brane solution will no longer be a solution with these varying param- 
eters. To be a solution, we have to add the corrections #mjv> an d X 1 ^ to the zero'th 
order solution with varying parameters, which should be chosen to satisfy the equations 
of motion. These corrections will be proportional to the first derivatives of the varying 
parameters, and we can neglect their variations along x M as it would be second order, so 
these corrections g^N> anc ^ X 1 ^ are functions only on the r-coordinate. One can 
choose the gauge using coordinate re-parametrization and gauge transformations to be 







9$ ~ % 



AW = 



(1) 







(3.18) 



Let us write the O'th order metric as 



ds 2 = -A(r)dt 2 + 2B(r)dtdr + C(r)(dx 



i\2 



(3.19) 



with 



A(r) = H~i(r)f(r) , B(r) = H~*(r) , C(r) = r 2 m(r) . (3.20) 
Then the metric up to first order in derivatives including the correction g^N looks as 



ds 2 



-A(r)dt 2 + 2B{r)dtdr + C(r)(dx 



i\2 



+ 
+ 2 



dt 2 + 2 



a? (dpB) + gil\r) dtdr + 2 [-x" (d^Ui) B(r)} drdx 1 



x» (d, Ul ) {A{r) - C{r)) + g$\r) 



dtdx 1 + 



x' l {d li ,C)S ij + g§ ) {r) dx i dx j (3.21 



and the gauge fields become 
Al _^mqj_ 



dt 



r* + qi 



+ 



^ ( + 4 (1, (r)l dt + U (d, Ui ) + 4 (1) (r) 

\r z + qi J [ r z + qi 



dx\{?>.22) 



4 The last gauge is different from the one in ref. [TH [17l HH [21] , but same as in ref. [HI [20] 
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and finally scalar fields will be 

x' = ^ + W^+*' (1) (0 • ( 3 - 23 ) 

Hi[r) v ' 

The task is to insert the above into the original equations of motion to obtain the equations 
for the first order corrections g$ N , A 1 ^ and X 1 ^, and then to solve them. 

Because we have a spatial 5*0(3) symmetry in our rest frame, one finds that the equa- 
tions for the first order corrections decompose into different 50(3) representations. After 
a tedious but straightforward calculation, one finds the following ordinary differential 
equations for g$ N , A 1 ^ and X 1 ^. The easiest one is the tensor mode, that is, traceless 
ij components of the Einstein equation, 

" Y r dr (r 3 f(r)d r (^^^j ) = Yr &r ^ H ^) + ^ ~ \ ^ S * 

(3.24) 

The vector mode equations are more complicated. From ti component of the Einstein 
equation, one has 

-oWf-A (r 5 H(r)O r ( -J?-) ) - £ %^ (O r A^) 
2r 3 H(r) \ V ; M r 2 #3(r) J J ^ r 3 H{r) \ ri ) 



f(r) ( 2m 



m(r) V 3 f(r) 
and from the ri-component, 



1 



2r 3 i^2(r) \ \r 2 Hi{r)J J i = \r z Hz{r) 

= + • (3-26) 

where is a complex expression in terms of first order spatial derivatives in g/s only, 
which is given in the Appendix. From the i-component of the Maxwell equation for each 
/, we have 

_ 1 ( r/(r)(ft(r))' { , m) \ _ ( _J^_ 

r \ H(r) V J J r y r 2 H^r) 

= i d , ( idtU ,)) + la r {-\c, JK ( {9j u k) ) 

r V rH*(r) / r \ 2 ( r + <u)\ r + Qk) J 

+ 1 -d r [ , \ (m(r 2 - qi ) (d iqi ) + qj(r 2 + Ql ) (d t m)) ) . (3.27) 
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The above five equations are vector mode equations. We mention that the ordinary 
differential operators in the left-hand side of the equations in the above and below take 
integrable forms, which is a quite non-trivial fact that has been checked via complicated 
algebra and educated guesses^. 

Finally, the most complicated part is the scalar mode equations under 5*0(3). From 
the tt-component of Einstein equation, one has 



-M) 9 - ( ff_i(r)/(r) ) * ( ff * (r)9 ' <;) ) - i /(r) y r) {h ^ )£ ' 

3r"ffi(r),.A r 2 + q, ) " 1 ' 



from the tr-component, 



^d r (r*Hhr)d r g$) + i V (d r A 



t 



2r 3 #5( r ) v y 3^r 3 i72(r) 

+iff*(r)* (^"'(r)/(r)) 9 r (i^r)^) + (#*(r)& ( 

E f (r 2 + g,) 2 + = Sg> , (3.29) 

fe! V r 2 + qi J 



4 /(r) 
3 r 4 #§M 



the rr-component looks as 
9 r (log (r>^(r))) a r (i^(r)<^) - £ r ( log ( ^ ) ) d r ( ) = 0(3.30) 



7=1 



and the last scalar mode equation from the Einstein equation is the trace part, that is 

E?=i(«), 

— d r (rm(r)d r (r 2 m(r)) g$) - - ]T . /^777 /t 7 « 
2r v r /=1 



—d r (rm(r)d r (r 2 m(r)) g$) - - £ v^7 (<Mt 



4 



r 2 # 



^7T^:(( r2+ "' 2 -^TV)• Y ' ,I, =^:s^ , . (s.sd 



5 In verifying these, we sometimes used Mathematica for basic algebra manipulations. 
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where and J2f=\ S[p are source terms proportional to space-time derivatives of 

black-brane parameters, which are given in the Appendix. 

One obtains more scalar mode equations from the Maxwell equations. From the t-part 
of the Maxwell equations for each /, one gets 

and from the r-component, 

r 3 H^(r) V J r 3 H^(r) V ' 

_ ±ymi d (HAri XI{1) \ 

r 3 H*(r) \H*(r) J 
Lastly, we present the scalar field equations. For this purpose, we choose X 1 and X 2 as 
independent variables with X 3 = x i x s ■ We have 

0, I - '//" , [...I llHr } \\glP 



(v/mgT (<9 r A t 1(1) ) - jmq z (d r A 3 t {1 




#i(r)# 2 2 (r)// ■ > r 2 Hs(r) 

d4r*f(rn(^X^ + J^X^)) 
1 K ' \H l s(r) 2Hl(r) J J 

r 2 | _ N | 2m( H \ ( H l( r ) X l(l)\ 

, . (r 2 + g 3 ) + tP%) ( %P-XW + ^X 2 ^) = S*>,(3.34) 
r 2 H*(r)Y (r 2 + q 3 ) 2 J \H*(r) H*(r) J 

and the similar equation with X 1 and X 2 interchanged. The source terms S 1 ^ and S 2 ^ 
can again be found in the Appendix. In the Appendix, we also sketch our method of 
computations for deriving these equations, the main task being to obtain the variation of 
Ricci tensor up to first order. 

The main point in this heavy endeavor is in fact to solve the above equations to find the 
first order transport coefficients; luckily enough, we are able to solve the above equations 
in explicit integral forms. 
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3.1 The solution 



We first observe that some combinations of the above equations are in fact constraints on 
the space-time derivatives of the black-brane parameters; they can't be arbitrary but have 
to be consistent with the conservations laws in the CFT side, such as energy-momentum 
and current conservations. In other words, bulk equations of motion include the conser- 
vation laws in the CFT side. Writing the Einstein equation as Emn an d the Maxwell 
equations as Mf f , one has three kinds of constraints equations, 



= g rt E tt + g rr E rt = -(HHr)S^ + Hs(r)f(r)Si 
{2mq 1 q 2 q 3 (<%«;) + md t (gig a g 3 )) 



J tr 



2r 9 Hl(r 



- ((d t m) (qiq 2 + qiqz + qm) ~ m dt (?i<?2 + <?2<?3 + ?3<?i)) r 2 

- (2m (diUi) (q x + q 2 + g 3 ) + 2 (d t m) (qi + q 2 + qs) - md t (q 1 + q 2 + q 3 )) r 4 

- (Am(diUi) + 3(d t m))r^ , (3.35) 

= g rt Ml + g"Ml = ~ 2 (d t (^m) + (d iUi ) ) , (3.36) 

r 6 H a (r) \ / 

= g rt E ti + g rr E ri = "* 0- {d.m) + 2m {d tUi )) . (3.37) 

r^H 3 (r) \^ J 

Note that for the above combinations of equations of motions, the radial differential 
operators cancel with each other to leave the above algebraic constraints. The constraints 
are uniquely solved by 

4 

{dtm) = --m(diUi) , (^m) = -4m(d t Ui) , 
2 

{d t qi) = --qi(diUi) or equivalently, d t (y/mqi) = -^Jmq 1 (d^) , (3.38) 

where the first two equations imply the zero 'th order energy momentum conservation, and 
the last equation is the conservation of U(l) 3 global symmetry currents. Indeed, applying 
AdS/CFT dictionary to our zero'th order black-brane solution, we have 

Jf ] = ^u^ Pl u» , (3.39) 

57TG5 

whose conservation laws in our rest-frame = (—1,0,0,0) are nothing but (13.381) . One 
generically obtains conservation laws of (k — l)'th order from the fc'th order equations of 
motion. 
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We next solve the remaining dynamical equations. It is easiest to solve the tensor 
mode equation (13.241) . Integrating it gives us 



9«V) 



r2Hi M(- 2 ^£ dr W + C «I dr '^7M +Ci ) ' (3 ' 40) 



with some constants Cy and C^, and we define 

&ij = - (diUj + djUi - ^ {d k u k ) S i:j \ . (3.41) 

One has to put C-j = as it is a non-normalizable mode. The CV, is uniquely determined 
to give a regular solution at the horizon r = rn where f(rn) = 0; note that the two 
integrals in the above logarithmically diverge near r = r#, and having a cancelation 
between the two for a finite result uniquely fixes 

C ij = 2r%Ht(r H )<T ij • (3.42) 

The vector mode equations ( 13.251) . ( 13.261) . and ( 13.271) are harder. As we already solved 
one constraint equation from (I3.25|) and (13.261) . we need to solve only (I3.26P and (13.271) . 
Let us first integrate (I3.27P once, which gives us 



H(r) V ' v \r 2 m(r) 
[d t Ui) + -CuK-^f- -f^- -e 3 (djU k ) (3.43) 



rH\{r) V 1 lJ 2 1J " (r 2 + qj)(r 2 + q K ) 
1 



m(r 2 - qj) (d iqi ) + gi (r 2 + q T ) (0 t m)) + C\ = Q« J(1) (r) + C\ 



with some integration constants C/, and then consider the horizon r = where f{rn) 
0. Imposing a regularity on and g*jp at r = r#, the first term drops and we have 



\={Q?\r H )+C!) . (3.44) 



" 2 H Hi{r H )] 2^qi 

The point is that the left-hand side is independent of J, so that C\ can not be arbitrary 
but has to take a form 

C! = -tf 1 \r H ) + -£^-C i , (3.45) 
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with only one degree of freedom d, which is nothing but the value of gu\ r H) at the 
horizon. The next step is to use the above f!3.43j) to replace (drA 1 ^ 1 ^ in the equation 
(13. 26f) to get a second order differential equation for gp only; 





- E " f 2 ( ^'V. K'M - Of 'fa) + Jjg-c, 

I=1 rf(r)(Hj(r)) \ r 2 H Hz(r H ) 

- r 2 ^( r ) L^) (^) + 2r 3 ^( r )p«( r ). 

To our surprise, the second order differential operator in the left-hand side is in fact 
integrable, that is, the left-hand side can be transformed into 

''WjAf^ll . (3.47) 




The way we have found this is the following. We start from the Ansatz for an integrable 
form, 

^d r (Qd r (S-)) = 1 [QSd 2 r ■ + (d r (QS) + Q (d r S)) d r • +<9 r (Q^S) ■) , (3.48) 
and comparing with the left-hand side of (I3.46p . one gets 

jQS = r 5 H(r) , i (0 P (QS) + Q (d r S)) = d r (r 5 #(r)) , (3.49) 
1 a Ma e\ V" 4mg 7 #(r) 

-d r {Qd r S) = 2 • (3.50) 

^ I=1 rf{r){Hi(r)) 

One can easily remove Q and S" in terms of P to get a differential equation for P, which 
turns out to be the same differential equation (I3.46p . but now without the source term 
in the right-hand side; that is, P is a homogeneous solution of the differential operator 
in (13.461) . Luckily, we know one way to generate a homogeneous solution without source 
terms; recall that any coordinate re-parametrization must correspond to a homogeneous 
solution of the problem. For our purpose, the following infinitesimal coordinate transfor- 
mation 

dt — > dt - edx 1 , dx % — ► dx l + e T dr , (3.51) 

r 2 H 2 (r) 



generates one homogeneous solution for which is H $(r)f(r), and from this one can 
let P be 

_ H-l(r)f(r) _ fjr) 

rmUr) rm(r) " ^ 



2 
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Once P is found, it is straightforward to find Q and S from the above equation to have 
the above integrable form of the differential operator. 

Solving (I3.46P then by integrating it once, we have 



aH ! («;M) V ii»«M 

M (£ i (9(Mi) + ^ V )) + c ; 

J oo 



H2 



with an integration constant C[, which can be fixed by considering the behavior at the 
horizon. Note that f(r) has an expansion near r = r# as 

f(r) = f'{r H ) (r - r H ) + \f"{r H ) (r-r H ) 2 + --- , (3.54) 

with f'{r H ) > 0, and the left-hand side in the above (13.531) takes a limit as r — > r#, 

r(f(r)f a ( Hl{r) r H f'(r H ) (i) r H f(r„) 

tt ( \ °r f ( \ <fc -» — ^r- — r9u ( r H) = — -17 — -Ci , 3.55 

where we have used the fact that g^\r H ) = Ci previously. This fixes C[ to be 

C t = -r*>IV) , (3.56) 



00 



where J(r') is the same integrand in (I3.53p . and one can rewrite f!3.53j) as 

«J^ 9 «U/V, M -^ C( . ,3.57) 
H{r) \ f(r) J Jr H Hz(r H ) 

We then integrate the above once more to have 

£\r) = (/W) - ^Mc\ + JgLc;, 

where one needs to put the integration constant C" to be zero as it corresponds precisely 
to the coordinate re-parametrization that we have used to get a homogeneous solution. 

At this point, the only uncertainty we have to fix is the constant Cj that appears both 
in J(r) and the above result. It might seem that it can be fixed for the to have a 

14 



regular derivative at the horizon r = r#; observe that the above result for g^' already has 
a finite value at the horizon irrespective of the constant Cj. Suppose that the integrand 
in the above equation has an expansion near r = r#, 

H ^ frV^.Is^^^^. +_»+... , (3.58) 



r' (/M) V™ -ff'M 7 (r'-r H ) 2 (r'-r ff ) 

then the first piece is harmless after integration as it cancels with f(r) in front, while the 
second piece will result in 

9u ~ ( r - t h) log(r - r#) H , (3.59) 

which has a divergent radial derivative that would signal divergent curvature tensors. 
Explicitly, one finds that the Ricci tensor R ri diverges with this. Therefore, we may have 
to choose right C, to make sure that b = in the near horizon expansion. Explicitly, one 
has 

„_ H(r H ) ( 1{rH) _M(^_ L _rp^ c \ t 



rH(f'(r„)f \' v "' Hi(r„) \H(r H ) r„ f'(r H 
and moreover one finds from (I3.53j) using f(rn) = that 

I(r H ) = V x ~ 4mgJ ,Q , (3.61) 

so that b is in fact proportional to C{ = gu\ r H), and it appears that we may have to 
put it zero for regularity. However, an explicit computation shows that the coefficient in 
front of Ci in fact vanishes identically, and the geometry is smooth for any value of C,. 
We need an extra input to fix the constant Cj. 

The answer to this puzzle lies in the frame choice, more concretely, the choice of either 
Landau frame or Eckart frame. In our work, we choose the Landau frame which states 
that 

Ujj T^ = -eu v , (3.62) 

and in particular, T t% = must hold in our local rest frame. In holographic renormalization 
that we will discuss in more detail in the next section, this condition gives us the constraint 
that the coefficient of \ in near boundary expansion of g^\r) should vanish, because it 
is precisely proportional to the first order correction to T t% . One can easily check that 
this indeed fixes our integration constant Cj uniquely. We have then completely solved 
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for Qu\r), as given above with Cj determined. Our explicit calculations give us 

4m V/tlVv^V 2r%m{r H )) J 



3 (^ + g/)(^ + gj)(r^ + g x ) / 

where is given in the next section. 

Once we have the solution for (r), one simply plugs it into the equation (I3.43|) to 
solve for A{ , whose integration gives us 

4« = r dr > - m. ( Q^y) - qf)( r .) _ 2^ ( oSV) c 



r'f(r') (#/(r')) 2 \ l 1 V \ (r')*H* (r') r 2 H Hl(r H ) 

(3.64) 

where we have chosen the integration constant to remove non-normalizable modes. This 
completes our long solution for the vector modes under 5*0(3). 

Finally we come to the solution of the scalar modes under 50(3). Although it seems 
difficult to our eyes to systematically solve these equations, we are lucky to be able to 
solve them by an educated guess from the previous analysis in ref. |16L 121] ; in fact, for 
the case of single U(l) R-charged hydrodynamics, most of the scalar modes under 50(3) 
turn out to be zero, except g# . Assuming the same feature, one easily finds that 

= %-H-^{r){d l u i ) , (3.65) 



indeed solves all the scalar mode equations of motion in the previous section_|. As the 
solution is expected to be unique up to trivial coordinate re-parametrizations, we conclude 
that the above g$ with 

gg> = 4 W = = , (3.66) 
is the solution of the scalar modes under 50(3). 

3.2 The first order transport coefficients 

It is a standard AdS/CFT procedure to obtain the first order corrections to the CFT 
energy-momentum tensor and the U(l) 3 symmetry currents from the results in the previ- 
ous section. We first discuss the energy-momentum tensor. One way to compute the CFT 



3 We used Mathematica for basic algebra manipulations in showing this. 
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energy-momentum tensor is to rewrite the full first-order metric in the Fefferman- Graham 
coordinateH where 

dp 2 

ds 2 = — — + p 2 giiv(p, x)dx' 1 dx u , (3.67) 
p 2 

and to read off the coefficient of the large p-expansion of g^ u (p, 



gi A )(x) 

g^u(p, x) ~ T]^ H + ^ +■■■ . (3.68) 



P 4 



The holographic renormalization procedure would then give us 



T ^ = ^G~/$W • (3-69) 



However, if one naively applies this to the zero'th order black-brane solution f)3.15p . one 
does not find the previously quoted energy- momentum tensor in (13.391) . 

T^® = — ^— (rf + 4uV) . (3.70) 

167TG5 

This is due to a subtlety in the scalar fields sector of the STU model; the scalar fields sector 
provides the cosmological constant at its vacuum, and hence the boundary counter-term 
that one adds in the holographic renormalization involves a non-trivial potential term of 
the scalar fields X 1 [38J. Because X 1 in the solution has a non-trivial profile, it turns out 
that this counter-term gives an additional contribution to the energy-momentum tensor 
other than (13. 69ft . The careful analysis including this subtlety was carried out in ref . |23j 
to find the above correct answer. 

For our purpose to find the first order correction T^ 1 ', we are however in a lucky sit- 
uation. Because the first order corrections to the scalar fields X 1 ^ vanish, there wouldn't 
be any first order contributions from the scalar fields sector to the energy-momentum ten- 
sor, and we can safely use (I3.69|) for the first order corrections to the energy-momentum 
tensor. Be warned that this may not be true in higher orders. A direct expansion of 
our first-order solutions in the Fefferman- Graham coordinate p which are related to our 
r variable by 



p + - + ^l^fiU +... , (3.71) 



( djUj) _ (Ejg/j (dm) (g j gj) 

3 6p 54p 2 

one finds that the only non-vanishing first-order correction to the energy-momentum 

comes from gQ' , and is given by 



Ti;> = ~2^f^U, = = -2 Wij , (3.72) 



7 It is also possible to get the energy-momentum tensor directly in the Eddington-Finkelstein coordi- 
nate, but the end results should be same in the first order in derivatives. 
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where we have used the fact that the horizon area per unit CFT volume is given by 
rjjH^(r}j), and the entropy density from the Bekenstein-Hawking formula gives 

S = ^g^l . (3.73) 

The last equality is simply the definition of shear viscosity rj, and one recovers the famous 
ratio 

- = h , (3-74) 

S 4:71 

in the STU model. It is not hard to make the previous result in a manifestly covariant 
standard form away from our static frame = (—1, 0, 0, 0); 

— p {rf v + Au tx u v ) — 2t7<7^ -\ , (3.75) 

with 

where P^ u = jf v + u^u u is the projection to the transverse components to 

It is also possible to write the first-order corrected metric in the covariant form, 

ds 2 = -H 2/3 {r)f{r) Ull u u dx»dx u - 2H 1/6 Ufl dx»dr + r 2 H 1/3 (r)P pu dx»dx v 
+ ^rH~ 1/6 (r)(d p u p )u fM u u dx^dx ,/ 

- 2 /(r) u ( f r dr' H{r>) ( [ r 'dr"lW(r") - r J^I^l C ) \dx»dx» 
H^(r) U "\L t'P{t') \L v 1 ] HW(r H ) *) J 

+ ^ (rj )^ W ^(£^(^-^^))^ .(,77) 
Here, the covariant first-order radial function ljp(r) is defined as 

+ 2 i$r) p >) • ( 3 - 78 > 

where the first-order functions Q^ir) and -P„^(r) are defined respectively as 
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The Cfj, is the covariantized form of C{ determined before. 

We next obtain the first order corrections to the U(l) 3 currents, which haven't been 
computed before in the literature, and would be the first non-trivial results in this work. 
The standard AdS/CFT formula 

J ? = }^l£^^Ai(p) , (3.80) 
works fine here, and it is easy to get the covariantized first order correction 

J F=vk{^™-^f-) ■ (3 - 81) 

where Qv^'ir) is defined in the above. Using the covariant version of the conservation 
(J33HD, 

P»d u m = -AmP^ [u a d a u u ) , (3.82) 

and the expression for the density pi in (13.391) . one can rewrite the result in a more 
suggestive form, up to first order in derivatives 

J? = Pl u» - VjP^D uPl + Cie up ^u u d pU<7 + ■■■ , (3.83) 

where the diffusion coefficient T>j is given by 

v {rl-qi) (3g4) 
2r* H H k 2(r H ) 

and the parity-violating coefficient (j (originated from the 5D Chern-Simons term) is 



32vrG 5 V 1Jn (rl + qj)(rl + qK ) 

3m JKL (r 2 H + qj)(r 2 H + q K )(r 2 H + q L ) J 

We finish this section by presenting the full covariant form of the gauge potential up to 
first order in derivatives; 

- 2^m^7 / dr / dr / dr I M (r - —j- — -C„ \dx^ , 

Joo [r'YHfir') Joo r"f\r")\Jr H Hs(r H ) J J 

where appropriate functions are defined previously. 
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4 Hydrodynamics with SU (2) in arbitrary dimensions 



Our next subject is to consider non-Abelian symmetry dynamics in hot hydrodynamic 
plasmas, motivated by the iso-spin SU(2)j dynamics in the QCD plasma. Although one 
can embed our bosonic action into the well-defined AdS±/ 'CFT 3 set-up of the Tri-Sasakian 
compactification of M-theory to AdS^, we perform our analysis in arbitrary dimensions 
envisioning that any system with non-Abelian symmetry would be described, at least 
approximately, by our model. In fact, the action we study has the simplest form one can 
imagine with gravity and gauge fields in AdS. However, in general dimensions other than 
n = 4 the connection to real QCD will no longer be our main motivation. 

We will consider (n + l)-dimensional gravity corresponding to n-dimensional CFT. 
We restrict ourselves to the cases of n > 3 only, as the n = 2 case seems peculiar in our 
results. Our action contains gravity with SU (2) gauge fields^ 

C = —±—(R + n(n-l)-FS rN F* MN ) , (4.87) 
167rG n+ i v 1 

with 

F a MN = d M A a N -d N A a M + e abc A b M A c N . (4.88) 

Note that we are allowed to choose and have chosen a specific normalization for the 
cosmological constant for simplicity, while the normalization of the gauge fields in the 
above corresponds to a definite value of coupling constant. In n = 3, this is dictated by 
the supersymmetry of N = 3 gauged supergravity [39], and we simply extend it to any 
dimensions. One can easily recover the gauge coupling constant dependence in our results 
below, if needed. The equations of motion are 

Rmn + (n + - — ^ Fp Q F aPQ ^j g MN — 2Fp M F aP N = , 

V M F aM N + e abc A b M F cM N = , (4.89) 

where Vm is the covariant derivative with metric Christofel connections. 
A charged black-brane solution in a general boosted frame is 

ds 2 = —r 2 V \r)u^u v dx^ dx v — 2u fM dx fl dr + r 2 (77^ + u^u u ) dx^dx v , 

A " = IS\7^ u > 4x * • (4 ' 90) 



8 One can simply substitute e abc below with any structure constants of a Lie algebra to have the general 
results for arbitrary Lie groups. 
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with 



V(r) = 1 - — + - 



.2n-2 



(4.91) 



where it is simply obtained by embedding the U(l) Reisner-Nordstrom black-brane into 
a Cartan direction inside SU(2) which is specified by q a (a = 1,2,3). As we are going 
to consider slow variations of q a over the CFT spacetime the Cartan U(l) will cor- 
respondingly vary point-by-point, and the non-Abelian nature will manifest itself when 
these variations are not parallel to q a locally. 

Let us consider slowly varying parameters u^, m, and q a up to first order, and we 
work in the frame where = (—1, 0, • • • , 0) at the position x M = 0. Then at first order 
in derivatives, we have 



and the above black-brane solution will no longer be a solution with these varying pa- 
rameters. To be a solution, we have to add corrections g$ N and A a ^P to the zero'th 
order solution with varying parameters, which should be chosen to satisfy the equations 
of motion. Our gauge choice is as before; 

n— 1 

9^ = , , 4^ = , £^ = . (4.93) 

i=i 

The resulting metric and the gauge fields at first order are 



u 



m 



m (0) _|_ x VQ^ m 
q a(0) + X „Q a 



(4.92) 



ds 2 



r 2 V^\r)dt 2 + 2dtdr + r 2 {dx 1 ) 2 



+ 




+ 



+ 



2 [a? {d^Ui) r 2 (v {0) (r) - l) + g£\rj\ dtdx 1 + 2 [-x" (d^m)} drdx 
g^Mr)dx l dx^ , 



(4.94) 



A 




(4.95) 
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After a lengthy calculation, we get the following equations for the first order corrections 
g$N an d ■ From the tt-part of Einstein equation, 

('•"" 1 *9. ( < 1) ) - \r'V(r)3 r (rV(r)) (a, 9 «) - 2n (r 2 V») 9 « 

.4 Mr v w)( ^-(^((M_^) , (4 , 6) 

from the tr-part, 

{r^drg^) + \d r (rV(r)) {d r g$) + 2n ft 



(i) 

tr 



from rr-part, 

^(S r9( «)=0 , (4.98) 



and from X)" =1 (ii)-part, 



^— T<9 r (r- 2 ^)) + (n - l)rV(r) (erf) + 2n(n - 1) 



2 (1) 

r 9tr 



-W (n " 1} 2 (n ~ 2) ^ (^f>) = 2(n - l)r (*«,) . (4.99) 



From the i-part of Maxwell equation, we have 



/(n-l)(n-2) 1 



2(n - 2) r 2 
r ((<9 t <f) + g a (<9^)) , (4.100) 



and from the r-part, 



The above six equations are scalar modes under SO(n — 1) spatial rotations. 

Vector modes equations of SO(n — 1) are the following. From the ti-part of Einstein 
equation, 

- m°> (#)) - s/^f^m m<") 

= l(M + ((^ rnr) + ^ (9tll ,) , (4 ,o 2 ) 
22 



from the ri-part, 



(#)) ) = , 

(4.103) 

and from the z-components of Maxwell equation, we have 



r n-3~' i r 2 i -V2fn-2)" — - r 1 1 



" ' ^ mq a ) + q a (d t u t )) . (4.104) 



V 2(n - 2) r 

Finally, the tensor mode, that is, traceless r/'-components of Einstein equation is 

(4.105) 

We now present the complete solution of the above equations. 



4.1 The solution 

We first solve the scalar mode equations. The constraint equations will be discussed after 
that. From (I4.98p . one finds that g^) = C with some constant C. To understand its 
meaning, note that this g^) = C will affect only (14.961) . (14.971) . (14.991) . and it is easy to 
check that it can be compensated by turning on 

g$ = -2Cr 2 , (4.106) 

that is, C corresponds to the above homogeneous solution of the problem. As the above 
g$ is a non-normalizable perturbation to the boundary CFT metric (look at the r 2 factor 
in front), we see that C in fact corresponds to a non-normalizable homogeneous solution 
of the problem, and we set it zero. Then equation (I4.10ip is integrated to give us 

Af ] = ^ + C a ' , (4.107) 

where C a> = is again a non-normalizable mode, and the meaning of C a can be easily 
understood by looking back the zero'th order profile of the gauge fields 

< (0) -i . (4-108) 
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9u = jz-^r (d lUi ) + — , (4.109) 



that is, C a is simply mapped to a redefinition of the charges q a } so that we can also set it 
zero. Then one can easily integrate (14.991) to have 

, — — r (diUi) + - 
[n — 1) r 

with an integration constant C . However, recall that the zero'th order g£' is 

(0) 2t r { \ 2 , 771 ( f ( f (a 1r >\ 

0W = _ r y( r ) = _ r j (4. 110 ) 

so that C is simply a redefinition of the energy density m. In summary, the only non- 
vanishing scalar mode in the solution is 

9$ = T -^-r S r(d i u i ) . (4.111) 
[n- 1) 

One then finds the following three kinds of constraints equations, 

= 9 rt E tt + g rr E r f = ((d t m) + j-^m (d^ - ((d t q«) + q« (d iUi )) 



= g H M? + g"'M: = sj {n ^ ((d t q a ) + <f (d iUi )) , 

= g rt E tl + g rr E n = --^ ri ((d l m)+nm(d t u t )) , (4.112) 

2r 

which results in 

(d t m) = -- m (Q. u .} (dim) = -nm (d t Ui) , (<9 t <f) = -g° . (4.113) 

(ra- 1) 

As before, they are the conservation laws for the zero'th order energy-momentum tensor 
and the SU (2) currents^ 

T M0) = ™ + nM ^) = p ^ + nufu 1 ') , 

167rG„, + i 



J^°) = —±—J ( n ^ 2) g VE P V . (4.114) 
47rG n+1 V 2 

We next solve for the vector modes, Eqs. fl4.103l) and (14. 104j) . which will turn out to 
have an important new ingredient due to the non-Abelian nature. An inspection shows 
that the only chance to see non-Abelian nature at this order is when (diq a ) is not parallel 
to q a . Define 



5 We will obtain them more rigorously in the next section. 
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then the following three SU(2) vectors in the Lie algebra of SU (2) form a normal basis; 

{q a ,Qle abc q b Q^} , (4.116) 
and especially one can expand (diq a ) in terms of them as 

(d,q a ) = (^W) q a + ? hc q h Q\ , (4.117) 

\ q-q J \ q-qj 

where P ■ Q = P a Q a . We then have to expand our in terms of them as 

Af 1] =A^\r)q a + f^(r)Q1 + g^(r)e abc q b Q c l , (4.118) 

with radial functions A^\ f( l \ and to be determined. Observe that we have dropped 
the index i for and as we will see that they are independent of ?o Plugging this 
expansion into our equations ( 14. 103f) and ( 14.1041) . one obtains the following two equations 
for Af> and gjp , 



(4.119) 



V 2(n - 2) r"- 1 \ <f • g 
and the following coupled equations for and g^ 1 ), 



- 1 1 tf" CM) + (flW ) , (4.120) 



^ (r»-V(r)a r /W) + 2] j 



Jl^l4dJC)-0, (4.121) 



2(n-2)r n_ s \ H 



1 A (^(OfW) - 1 9r f Zf) = 11 (4 .122) 



r"- 3 v w ' f2(n-2)rH \r5 / V 2 (™ ~ 2 ) r " _1 ?' 9 

It is straightforward to solve (14. 1 19[) and (14.1201) as is done in the STU model, and we 
simply present the result; 

ft'," = rV(r) fdr' -— 1— , ( f dr" J(r") - rg- V(r,)C.) , (4.123) 



r y w JL ™ U. *'" /(r " ) " rr y MC ' 



10 This should be clear from the spatial 50(3) symmetry because Q1 and e abc q h Q1 are already SO(3) 
vectors. 
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where 



n ( \ I n-1 1 / q • ($g) 

Qt(r) = -W^ tt- ^ +(d t Ui) 

y 2(n — 2J r V g • q 



(n - l )(w-2) g- q ( „ ^ ^ , , , /(n - l)(n - 2) Ci 
2 r' 



J(r) = -4,/ ^ — jdl I QiW - Qi(rg) + 



(n-l)r"-V(r) (^) . 



Again, the integration constant Cf is fixed by the Landau frame constraint, and an explicit 
computation gives us 

= — "7 ^=T ^-^+ [ nr H „x n-l ^ • 

nm \ (n — 2)r H \ (n — 2)r# / / 

However, we are unable to integrate the equations (I4.12ip and (I4.122p to solve and 
g^ 1 '; we will instead comment on a possible numerical approach. Defining f = 
+ i\q\g^ l \ the equations (I4.12ip and (14.1221) become a single complex equation 

-\d r (r^V(r)d r F^) - 2i\g\J J ~ \ 'd r ( = -i\J ^~\ J-J-_ 

The equation has one trivial solution, 

F (+)(D = L (4.124) 

9-9 

which is non-normalizable. However, having this solution is of great help in finding 
the unique normalizable solution regular at the horizon numerically; one only solves the 
homogeneous equation without the source term in the right-hand side, and then add 
the above trivial solution to have a normalizable solution. Considering the limit of the 
homogeneous equation to the horizon where V(th) = 0, one gets the relation 



{d r F^ 



"^l^ 1 , (4.125) 



r =r H ^2(n - 2)r£+V(r H ) - 2i^/n~=l\q\i 

for a regular homogeneous solution. Putting [r h) = 1, the above completely speci- 

fies the boundary condition at the horizon, and one can numerically solve the differential 
equation uniquely, that is, regular homogeneous solution normalized asFWW(r H ) = 1 is 
unique. Let's call this homogeneous solution Fq + ^(t). In general, its large r asymptotic 
will give us a finite non-zero constant ^(oo) 7^ 0, and one constructs the full solution 
of our original problem simply as 

^-Tfi'-W^l ' (4 ' 126) 
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Finally, the tensor mode equation (I4.105P is easily integrated to give us 

g$> = -2,„r" £ dr> (r<) „ + 1 1|/(r , ) (M- - (nr)- 1 ) , (4.127) 

with 

= ^ (W,) + - . (4.128) 

4.2 The first order transport coefficients 

Based on the results of the previous section, we can write the covariant form of the metric 
and the gauge field up to first order in derivative expansion. The metric looks as 

2 

ds 2 = —r 2 V{r)u p u y dx^'dx u — 2u p dx iJ 'dr + r 2 P pv dx^dx u -\ -r(d p u p )u p u u dx^dx 1 ' 

~ r2V ^{lj r \ r >)n-\ V 2 { r>) [Syr''W')-r^V'{r H )c)i )dx»dx» 

/ rr 1 \ 

„2_ / / j l I IJ\n—\ /„ \n-l \ \ j„u j„v 



2r a ^ { L dr> {r')^V{r') W ~ {TH) ) ) dx " dxV ' (4 ' 129) 
where the first-order functions I^p{r) and Q^\r) are defined respectively as 



+ (n - l)r n ~V(r)M^M 



and we also defined 



a ,u = _ pmp „(3 {QaUp + dpUa) _ _____ {daUa) (4 131) 



The gauge field is found as 



/(« - i)( B - 2) q,\ _ /(«-!)(« -2) /V 1 



2 r?J V 2 7oo fr')"" 1 



r." 



dr \ rr+ l v2{rM) (£ *'"/f>n-4-V'(r«)^| V 



+ (/ (1) (r)Q«+^ 1 )(r)e ofc g b Q«)d^ 
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where = e q P"{d v q c ). Note that the last line is the non-Abelian induced terms. 

Again using holographic renormalization, one can easily find the stress tensor and the 
SU(2) charge current. In terms of Fefferman-Graham coordinate expansion, 



9M = ^u + --- + ^ + --- , (4.132) 



the n-dimensional CFT stress tensor is given by 



r ~ = idb*' • (4 - 133) 



Again, it is straightforward to check that the only non-vanishing first order correction 
comes from gff with 

T^ = -2 V a i3 , (4.134) 

where 

r n-l „ 
lo7rG n+ i 47T 

We point out that it is not a trivial fact here and in the STU model that there is no first 
order correction to T 00 , as it results from a non-trivial cancellations in the expansion. 
The full covariant expression of the stress tensor is 

T» v = p (jf + nu»u v ) - Tiqa^ + •• • . (4.136) 

with the pressure p = 167r ^ representing the zero'th order contribution. 
The SU(2) charge current is similarly obtained from the expansion as 

J >hd^ pn ~ 2AW • (4 - 137) 

and it is easy to find the first order correction to be 



where (r) is given before in fT4TT30|) . and f {n ~ 2 \ 

g(n 2) are coe fg c i en ts of l/p n 2 in 
the expansion of f^'ip) and g^(p) respectively. Trading q a with the density p a defined 
in (14.1141) . we can rewrite the result in the form 

jail) = _ v ^ • Jfrp) _ u v djjU \ p a + p i£ ^( 9/ ) 

+ v 2 p;{p a {p-d uP )-{ P - P ){d uP a )) , (4.138) 
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with three diffusion coefficients 

v = 1 L _ 2q-q \ = (n - 2)m + 2r n H 

(n — 2)r H \ nmr 1 ^ 2 J n(n — 2)mr H 

Vl ~ (n-1) ' P2 -(n-l)l(n-2)l ■ (4 ' 139) 
The 22 is essentially the usual diffusion coefficient of Abelian nature, which agrees with 
ref.jlOj, while the other two diffusion coefficients are due to the non- Abelian properties. 
Although their precise values can only be determined by numerical analysis, we hope that 
the above structure of non-Abelian current we obtain in derivative expansion may be im- 
portant for future applications to the QCD plasma incorporating non-Abelian symmetry. 



4.3 On Tri-Sasakian compact ificat ion of M-theory to AdS^ 

We would like to conclude by a few comments on the realization of our SU(2) theory 
in a concrete example. Our 4-dimensional bosonic action (n = 3) with SU(2) gauge 
symmetry in the bulk has a specific AdS^/CFT^ realization in M-theory; consider N 
M2 branes sitting at the apex of an 8-dimensional Hyper-Kahler cone and take a near 
horizon limit. The superconformal theory one gets on the M2 branes has iV = 3 or 6-real 
components of dynamical supersymmetry with SU(2)r R-symmetry. The corresponding 
dual theory on AdS± will include a consistent truncation to the minimal N = 3 577(2) 
gauged supergravity in 4- dimensions [39], whose bosonic action is precisely our action in 
the previous section. In this case, one has an explicit expression for the Newton's constant 
G4 in terms of the number of M2 branes as follows. 

The 11-dimensional M-theory supergravity action is 

c ^^J dUx ^( R{n, -\^ 2 )-idmJ c ^ F ^ F ' • (4 - 140) 

where is the Ricci scalar of the 11-dimensional metric, and 

\F 4 \ 2 = ±F MNPQ F MNp Q . (4.141) 

The near horizon limit of M2 branes at the tip of a Hyper-Kahler cone takes a form of 
AdS± x X7 with X-j being a Tri-Sasakian 7-fold which is the unit radius section of the 
Hyper-Kahler cone involved. The explicit solution is given as 

dsn = R 2 {\ds 2 AdS4 + dn 2 x ^j 
F 4 = ^i? 3 6 4 , (4.142) 
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where dVL 2 x is the metric of X-j normalized in such a way that 

Rab = Q9ab , (4.143) 

and €4 is the volume form of the unit radius AdS±. The constant R is given by the relation 

QR 6 vo\(X 7 ) = {2itl p fN . (4.144) 

Then one can easily obtain the 4-dimensional effective action on AdS± after compactifying 
M-theory action (14.1401) on the above 7-dimensional Tri-Sasakian manifold X7. Note that 
the X-j metric of R 2 dQ 2 Kr now has 

Ra b 7 = ^9aJ , (4-145) 



so that R x? = jp, which means 



Also one has 



R (n)^ R (*) + ^ . ( 4.i4 6) 
Br 



^ 2 (4V=S ■ (4-147) 



\R 2 J R 2 



Combined with 



d n x = i? 7 vol(X 7 ) J d 4 x y/=gl , (4.148) 

the effective 4-dimensional action takes a form 

T^f /<W=*(*» + | + -) ■ <«*> 

and we identify 



i? 7 vol(X 
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1QtiG 4 (2vr) 8 /9 



(4.150) 



We then need toput R = 2 to conform to our convention of cosmological constant in the 
previous section_J, and using (14.1441) one finally has 

1 2vriV§ 

(4.151) 



167rG 4 2 2 6§ (vol(X 7 )) 5 

For a class of Tri-Sasakian manifolds that are obtained from Hyper-Kahler quotients, 
their normalized volumes vol(Xj) are explicitly known [4"Tj |4"2"] . One can start from a 



n The quicker way to arrive at this conclusion is to make the AdS± in (|4.142|) to have unit radius, which 
should be a solution with our cosmological constant convention. 
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(2 + r)-dimensional flat quaternion space and take U(l) r Hyper-Kahler quotients specified 
by charges Q l a where i runs over U(l) and a runs over (2 + r) quaternions. The resulting 8- 
dimensional Hyper-Kahler cone will have a unit radius section as a Tri-Sasakian manifold, 
whose normalized volume is known by the formula [12] 

or+l 7r 4 „ r 2+r i 
WOl(X 7 ) = — — —— r r / TTd^ IT o ■ (4-152) 

1 7> r(4)Voi (U(iy) J l\ v u i + (E[ =1 QL0 l ) 

In the simplest example of r = 1 with three charges Qi (i = 1, 2, 3), it becomes 

i fv (r\ r\ r\ \\ ^ (Q1Q2 + Q2Q3 + Q3Q1) (A , kq n 

vol (Xr(Q h Q 2 , Q,)) = Y (Q 1 4. Q 2 ) (Q 2 + Q 3 )(Q 3 + 0l) ■ < 4 ' 153 ) 

which includes the famous iV(l, 1) as Qi = Q2 = Q3 = 1 with 

vol(JV(l, 1)) = — . (4.154) 



Although the N = 3 superconformal theory on M2 branes dual to the gravity background 
with N(l, 1) is still unknown, there is a proposal in ref. [43] for the case of Qi±Q2±Q3 = 0, 
inspired by the theory of BL/ABJM|44[ 145] . See also ref.[46] for N = 2 cases. 
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A Appendix 1 

A convenient expression for the first order variation of the Ricci tensor we use is 

SRmn = -VmC^ + V P Cl IN , (A.155) 

with 

c mn = 7>9 PQ (Vm^jvq + V N 8g MQ - V Q 5g MN ) , (A. 156) 

where Vm is the covariant derivative with respect to the zero'th order metric gMN- The 
tensor C^ N is in fact a variation of the metric Christoffel symbol, C[ 4N = 6T^, IN . 
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B Appendix 2 



We denote 

A(r) = H-*(r)f(r) , 5(r) = fH(r) , C(r)=r 2 H^(r) , (B.157) 

as before, where the local parameters m and qi are implicit in the expressions above and 
below. The prime in the equations below means the radial derivative 4-. 



_1 Aft* 1 ) l (d t X*) \ _ ( 1_ ( (dtX 1 ) l (d t X 
B \ X 1 2 X 2 J \B { X 1 2 X 2 
1 (B' 30} ( (d t Xi) l (d t X 2 ) } 3 ((X 1 )' 1(X 



(diUi 



2\i- 



B \ B 2C J \ X 1 2 X 2 J 2BC v ' \ X 1 2 X 2 
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